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ABSTRACT 

A steady two-dimensional MHD free convection and mass transfer flow past an inclined vertical surface in the 
presence of heat generation and a porous medium have been studied numerically when the fluid viscosity and 
thermal conductivity are assumed to be vary as inverse linear function of temperature. The governing partial 
differential equations are reduced to a system of ordinary differential equations by introducing similarity 
transformations. The non-linear similarity equations are solved numerically by applying the Runge-Kutta 
method of fourth order with shooting technique. The numerical results are presented graphically to illustrate 
influence of different values of the parameters on the velocity, temperature and concentration profiles. Skin 
friction, Nusselt number and Sherwood number are also completed and presented in tabular form. 

Keywords - Heat generation, inclined vertical surface, mass transfer, shooting method, variable viscosity and 
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I. INTRODUCTION 

Under the influence of magnetic field, the study of free convection and mass transfer flow of an electrically 
conducting fluid over an inclined vertical surface has considerable interest in geophysics, astrophysics and many 
engineering problems, such as cooling of nuclear reactors, the boundary layer control in aerodynamics. Alam 
et.al [1] studied the heat and mass transfer in MHD free convection flow over an inclined plate with hall current. 
Anghel et.al [2] investigated a numerical solution of free convection flow past an inclined surface. Subhakar et. 
al [3] studied the Soret and Dufour effects on MHD free convection heat and mass transfer flow over a 
stretching vertical plate with suction and heat source/sink. 

The study of the flow through porous medium has some important applications in engineering and geophysics 
such as chemical engineering for filtration and purifications process, agriculture engineering to study the 
underground water resources, and petroleum technology to study the movement of natural gas, oil and water 
through the oil reservoirs. Chandra et. al [4] studied heat and mass transfer along an isothermal vertical porous 
plate in the presence of heat sink. Sharma et. al. [5] studied the effects of dusty viscous fluid on MHD free 
convection flow with heat and mass transfer past a vertical porous plate 

The study of the heat generation or absorption in moving fluids is important problems dealing with 
chemical reactions and those concerned with dissociating fluids. Reddy [6] studied the effects of heat generation 
and radiation on steady MHD free convective flow of micro polar fluid past a moving surface and they [7] also 
observed the MHD free convection heat and mass transfer flow past an inclined vertical surface in a porous 
medium with heat generation. Islam et. al [8] discussed the MHD free convection and mass transfer flow with 
heat generation through an inclined plate. 

Most of these studies are based on constant physical properties. More accurate prediction for the flow and 
heat transfer can be determined by variation of these properties with temperature. Several investigators have 
studied the effects of temperature dependent viscosity and thermal conductivity properties of fluid. But 
temperature dependent physical properties like viscosity of the fluid and thermal conductivity plays a significant 
role in fluid mechanics. It plays a great role in underground storage system and geothermal energy extraction. 
Khound et. al. [9] observed that a significant variation of velocity distribution and temperature distribution take 
place with the variation of viscosity parameter and thermal conductivity parameter. 
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Mahanti et.al. [10] studied the effects of varying viscosity and thermal conductivity on steady free 
convective flow and heat transfer along an isothermal vertical plate in the presence of heat sink. They assumed 
the viscosity and thermal conductivity to vary linearly with temperature. Sarma et. al [11] studied the effect of 
variable viscosity and thermal conductivity on heat and mass transfer flow along a vertical plate in the presence 
magnetic field and they also [12] studied on free convective heat and mass transfer flow with constant heat flux 
through a porous medium. 

The purpose of this present paper is to investigate the effects of variable viscosity and thermal 
conductivity on MHD free convection and mass transfer flow past an inclined vertical plate in a porous medium 
with heat generation. The flow governing equations are transformed into ordinary differential equations, before 
numerically solve by Runge-Kutta shooting method by using similarity transformation. The fluid viscosity and 
thermal conductivity are taken as inverse linear functions of temperature. The effects of variable viscosity 
parameter, variable thermal conductivity parameter, Prandtl number, magnetic parameter and Schmidt number 
have been discussed and shown graphically. 


II. MATHEMATICAL FORMULATION 

Consider a two dimensional steady laminar boundary layer flow of an incompressible viscous electrically 
conducting fluid over a semi-infinite inclined vertical plate in a porous medium with an acute angle a . The flow 
is assumed to be in the v -direction, which is taken along the semi-infinite inclined porous plate and y - axis 
normal to it. A uniform magnetic field B 0 is imposed along the y-axis. The fluid properties are assumed to be 


isotropic and constant except for the fluid viscosity and thermal conductivity. 

The flow governing equations for the present problem are: 

Equation of continuity 

du dv 

+ — = 0 

dx dy 


Equation of momentum 

du du 1 d 

u 1- v — = 

dx dy dy 


r du ^ 

u — 

dy. 


+ gP* ( T - T oo )cos(a) + gP** ( C-CJ cos (a) 
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Equation of energy 

f dT dT 

pc n u l-v — 
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+ QPT-TJ 
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(3) 


(Equation of concentration 

dC dC „ d 2 C 

dx dry dy 


(4) 


where (u, v) is the fluid velocity, //, p x ,k are the coefficient of viscosity, density and thermal conductivity of 
fluid, g , p and p are the acceleration due to gravity, thermal expansion coefficient and concentration 
expansion coefficient. k x , T and are the permeability of porous medium, fluid temperature within the 
boundary layer and fluid temperature in the free stream. C and D are the concentration and the coefficient of 
mass diffusivity of the fluid. Q x , C p are the heat generation constant and specific heat at constant pressure. 

The boundary conditions for the problem are: 

u = 0,v = 0,T = T W , C = C w at y = 0, 

u —>0,T —> T rj , C p C rt at y ^ go 

Lai and Kulacki [13] has assumed the fluid viscosity as 
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or, 
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E E r 


= c(T - T k ) 
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U ^ rp r 1 

here, c = —,l k = l r - — 

K % 

where // OT , k f and I , are the viscosity, thermal conductivity and temperature at free stream, a, c, Tr and T k 
are constants and their values depend on the reference state and thermal property of the fluid. In general a > 0 
for liquids and a < 0 for gases, y and f are constant based on thermal property of the fluid. 

Let us introduce the following transformations: 


U = UJ'(ri), V = \ A p^(/zf -f),rj= P -y , 
2 V x \ v x x 


( 8 ) 


T -T 

W GO 


C -C 

W 00 


Using the transformations (5-8), in equations (1-4), it is seen that the equation of continuity satisfies 
identically and rest of the equations becomes: 

or 
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9" - 


— — — i — ff" + GrO cos a + Gm(j) cos a — {M + K)f 1 = 0 
9 — 9 „ 9 „ 2 
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= 0 
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H 
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, <jB0x 

M = is the magnetic field parameter. 

pjo* 

Gr = ^-l— j s (p, e i oca i Grashof number. 


U, 


Gm 




u. 


is the mass Grashof number. 


jUC p 

P, = is the Prandtl number. 
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Sc = is the Schmidt number of fluid. 


0 


Tr-T* 

T-T 


1 


T -T 

Q _ L _k J Q0 


r(T„-Tj 

1 

~WPC) 


is the viscosity variation parameter. 


is the thermal conductivity variation parameter. 


K = 


k,U x 


is the permeability parameter. 


Q x x 

Q = is the heat generation parameter. 


PooCpU 


p oo 


Re = 


Ux 


v„ 


is the local Reynolds number. 


( 12 ) 


The transform boundary conditions are 

/ = 0, /' = 0, 0 = 1 , (ft = 1 at 77 = 0 1 

f' = 0, 0 = 0, 0 = 0 as r/ — >oo j 

The physical quantities of interest are the skin friction coefficient Cy , Nusselt number Nil and Sherwood 
number Sh for the surface which are define by 

0. ~ ~ ^ 




Nu = 


Sh = 


P,U y 


(e-e r ) 


Re 3 /'( 0) 


XC l, 
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xh 


D(C-C„) 
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where T - (J, 


du N 

.^Jv=l 


,P w = k 


dT ^ 


y=0 


\dy Jy=Q 


, h =-D 


v^y y=0 


III. RESULTS AND DISCUSSION 

The boundary value problems (9)- (13) are solved numerically using fourth order Runge-Kutta shooting 
methods. The numerical values of different parameters are taken as Pr = 0.7, d r -A,G k = 4 , M = 0.5, 
Gr = 0.5, Gm — 0.5, CC — 30 unless otherwise stated. Both velocity, temperature and concentration profiles 
are displayed graphically. 

Velocity profile for various combination of parameters Pr=.7, Ec=.05, N = 0.5 ,k = 0.2 are illustrated 
in figure 1 to figure 8 . It is observed from these figures that the velocity distribution decreases with increase of 
viscosity variation parameter 6 r , thermal conductivity parameter 6 k , magnetic parameter M ? permeability 

parameter K and Schmidt number Sc . On the other hand velocity increases with increasing values of heat 
generation parameter Q , local Grashof number Gr and mass Grashof number Gm. 

From figure 9 to figure 16, it can be concluded that the temperature increases with the increase of viscosity 
variation parameter 9 r , magnetic field parameter M , heat generation parameter Q and Schmidt number Sc. 
Also temperature decreases with the increase of thermal conductivity parameter 6 k , local Grashof number Gr 
and mass Grashof number Gm. 
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Concentration profiles for various values of Sc and M are shown in figure 17 to 18. It is clear from these 
figures that the concentration decreases with increase of Schmidt number Sc. But it is increases with increase of 
magnetic parameter M . 

The missing values f"( 0), 6^(o), ^(())and the coefficient of skin friction Cf , Nusselt number 

Nu and Sherwood number Sh are estimated for various combinations of parameters and presented in Table 1 
to Table 3. From these Tables we have observed that Cf and Nusselt number Nu decrease and Sh increases 

for the increasing values of viscosity parameter 0 r while f"{0) and Sh increase and Nu decreases with the 
increase of thermal conductivity parameter G k . Also Sh increases with increase mass Grashof number Gm and 
decreases with Schmidt number Sc . 

In the Table 4 we have compared of missing values /”( 0), 6^(0) and ^(o) with those of values obtained by 
Reddy et. al [7] for various values of Gr . From this Table we have observed that the accuracy of missing values 
/"(0), #'(0) and ^'(o) of present study are less than that of Reddy et. al [7]. 

IV. Figures and Tables 




Fig 1 : Velocity profile for different 0 r . Fig 2: Velocity profile different 0 k . 




Fig 3: Velocity profile for different A/. Fig 4: Velocity profile for different Q. 




Fig 5: Velocity profile for different Gr. Fig 6: Velocity profile for different Gm. 
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Fig 7: "Velocity profile for different •S'c. Fig 8: Velocity profile for different AT. 




Fig 9: Temperature profile for different 0 r . Fig 10: Temperature profile for different 0 k . 




Fig 12: Temperature profile for different M . 




Fig 14: Temperature profile for different Gm. 




Fig 15: Temperature profile for different Sc. Fig 16: Temperature profile for different K. 
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Fig 17: Concentration profile for different Sc . Fig 18: Concentration profile for different A/. 


Table: 1 


M 

0 r 

/'( o) 

0lo) 

f(0) 

Cf 

Nu 

Sh 


2 

3.862029 

-0.45512 

-0.22722 

438.6968 

-0.14488 

0.071854 

0.25 

4 

3.595297 

-0.4799 

-0.24603 

181.3049 

-0.15277 

0.077801 


6 

2.713935 

-0.48647 

-0.25091 

34.10165 

-0.15485 

0.079345 


2 

3.653211 

-0.44259 

-0.21555 

415.1085 

-0.14076 

0.068162 

0.5 

4 

3.405082 

-0.46589 

-0.23333 

171.7936 

-0.14817 

0.073787 


6 

2.582761 

-0.47205 

-0.23795 

32.48405 

-0.15012 

0.075245 


2 

3.469272 

-0.43109 

-0.20469 

394.3685 

-0.13694 

0.064728 

0.75 

4 

3.237415 

-0.453 

-0.2215 

163.4325 

-0.1439 

0.070043 


6 

2.466796 

-0.45879 

-0.22584 

31.06169 

-0.14572 

0.071418 


Table 2: Table: 2 


Q 

o k 

r(o) 

si o) 

/(o) 

Cf 

Nu 

Sh 


2 

2.355011 

-0.43614 

-0.25368 

29.63691 

-0.13337 

0.08022 

0.1 

4 

2.388106 

-0.42422 

-0.32137 

30.036 

-0.13453 

0.101625 


6 

2.513231 

-0.42103 

-0.3409 

31.51908 

-0.13911 

0.107802 


2 

2.432502 

-0.44259 

-0.21555 

30.67494 

-0.13548 

0.068162 

0.2 

4 

2.463873 

-0.43107 

-0.2688 

31.05703 

-0.13656 

0.085002 


6 

2.582761 

-0.42799 

-0.2838 

32.48405 

-0.14076 

0.089746 


2 

2.518574 

-0.44936 

-0.17475 

31.8122 

-0.13774 

0.05526 

0.3 

4 

2.547707 

-0.43832 

-0.21228 

32.16999 

-0.13873 

0.067129 


6 

2.658662 

-0.43536 

-0.22232 

33.51574 

-0.14258 

0.070302 


Table: 3 


Gm 

Sc 

/'( 0) 

si o) 

tho) 

Cf 

Nu 

Sh 


0.2 

2.514299 

-0.28718 

-0.21369 

31.624 

-0.09133 

0.067576 

0.25 

0.4 

2.506848 

-0.36825 

-0.21264 

31.53248 

-0.11711 

0.067241 


0.6 

2.500946 

-0.43836 

-0.21187 

31.45971 

-0.13939 

0.067001 


0.2 

2.609755 

-0.28969 

-0.21915 

32.81628 

-0.09215 

0.069302 

0.5 

0.4 

2.594611 

-0.37199 

-0.21705 

32.63021 

-0.11832 

0.068637 


0.6 

2.582761 

-0.44259 

-0.21555 

32.48405 

-0.14076 

0.068162 


0.2 

2.704833 

-0.29217 

-0.22448 

34.00385 

-0.09296 

0.070987 

0.75 

0.4 

2.681772 

-0.37565 

-0.22135 

33.7204 

-0.1195 

0.069996 


0.6 

2.663932 

-0.44671 

-0.21913 

33.50032 

-0.14209 

0.069293 
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Table: 4 


Gr 

/'( 0) 

0'(o) 

f(0) 

/'( o) 

ff(d) 



Ready et.al [7] 


Present study 


1 

1.355485 

-0.31921 

-0.12004 

1.048356 

-0.21489 

-0.07232 

2 

2.058302 

-0.36851 

-0.21732 

1.64646 

-0.23438 

-0.12823 

3 

2.695334 

-0.40412 

-0.27765 

2.192476 

-0.2491 

-0.16583 


V. Conclusion 

From the above analysis, it is clear that the viscosity and thermal conductivity variation parameter, heat 
generation parameter Q , magnetic parameter M , permeability parameter^, local Grashof number Gr, mass 

Grashof number Gm and Schmidt number etc. have substantial effects on velocity, temperature and 
concentration profiles within the boundary layer. We can conclude from this analysis as: 

1 . The increasing values of viscosity retard the velocity but enhance the temperature. 

2. Both velocity and temperature increases as heat generation parameter increases. 

3. Velocity decreases and temperature increases with increase magnetic parameter. 

4. Velocity increases and temperature decreases with increase of free convection parameter. 

5. Skin friction and Nusselt number decrease with increase of viscosity. 

6. Increasing values of Schmidt number retards the Sherwood number. 
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